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Abstract 

Saproxylic insects play an important part in decomposing dead wood in healthy forest ecosystems, but little is known about 
their role in the aftermath of large-scale forest mortality caused by pest insect outbreaks. We used window traps to study 
short-term changes in the abundance and community structure of saproxylic beetles following extensive mortality of 
mountain birch in sub-arctic northern Norway caused by an outbreak of geometrid moths. Three to five years after the 
outbreak, the proportion of obligate saproxylic individuals in the beetle community was roughly 10% higher in forest 
damaged by the outbreak than in undamaged forest. This was mainly due to two early-successional saproxylic beetle 
species. Facultative saproxylic beetles showed no consistent differences between damaged and undamaged forest. These 
findings would suggest a weak numerical response of the saproxylic beetle community to the dead wood left by the 
outbreak. We suggest that species-specific preferences for certain wood decay stages may limit the number of saproxylic 
species that respond numerically to an outbreak at a particular time, and that increases in responding species may be 
constrained by limitations to the amount of dead wood that can be exploited within a given timeframe (i.e. satiation 
effects). Low diversity of beetle species or slow development of larvae in our cold sub-arctic study region may also limit 
numerical responses. Our study suggests that saproxylic beetles, owing to weak numerical responses, may so far have 
played a minor role in decomposing the vast quantities of dead wood left by the moth outbreak. 
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Introduction 

Insect outbreaks are major disturbance factors in many forest 
ecosystems, periodically causing defoliation and mortality of trees 
and other vegetation across vast areas [1]. The potential causes of 
outbreaks have been intensively studied for decades. More 
recently, however, there has also been increasing interest in their 
ecological consequences [2]. One important consequence of 
outbreaks is that they alter resource availability for many 
consumer species. For instance, outbreaks can increase food 
availability for insectivorous birds [3-5], small mammals [6,7] and 
invertebrates [8,9]. Similarly, understory vegetation may benefit 
from increased nutrient inputs in the form of frass and insect 
cadavers [10], and increased availability of light due to removal of 
foliage in the canopy [11,12]. The resulting increase in understory 
plant biomass may, in turn, provide resources for species that 
depend on understory vegetation for food or cover [13,14]. In all 
of these cases, outbreaks provide consumers with resource pulses, 
i.e. rare and ephemeral events of resource superabundance 
[15,16]. For most systems, we still have very limited understanding 



of how such outbreak-induced resource pulses affect the dynamics 
of consumer communities. 

Saproxylic insects are one group of organisms that may 
experience particularly dramatic increases in resource availability 
following outbreaks. Saproxylic insects are, per definition, 
"dependent, during some part of their life cycle, upon the dead 
or dying wood of moribund or dead trees, or upon wood- 
inhabiting fungi or the presence of other saproxylics" [17]. Thus, 
since outbreaks kiU or damage a large number of trees within a 
short time, they may provide saproxylic insects with resource 
pulses in the form of massive amounts of dead or dying woody 
material. Accordingly, saproxylic insects could be expected to 
increase rapidly in abundance after outbreaks. Indeed, it is well 
documented that saproxylic insects respond numerically to large- 
scale forest mortality due to other disturbances, like fire [18-20] 
and wind [21-23], and recently this has also been reported after 
outbreaks by tree-killing bark beetie in central Europe [24]. 
However, the magnitude of such responses is not necessarily a 
simple function of resource input. One important limiting factor 
may be the intimate relationship between saproxylic insects and 
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the stage of decay of dead wood. A given saproxylic species can 
usually exploit dead wood only in a specific stage of decay. Thus, 
at any point in time after an outbreak, only a subset of the 
saproxylic community will be able to utilize the dead wood 
resource pulse. Mcjreover, the decay stage a saproxylic species can 
utilize is often relatively short in duration, lasting no more than a 
few years [25-27]. This is likely to lead to satiation effects, where a 
species is able to colonize and exploit only a small fraction of the 
dead wood generated by an outbreak within the time period when 
the wood represents a usable resource. Consequendy, numerical 
responses to outbreaks both in individual species and the 
saproxylic community as a whole may in many cases be relatively 
limited in magnitude, compared to the magnitude of dead wood 
input. It should also be kept in mind that insect populations may 
be naturally limited or regulated by numerous other abiotic and 
biotic factors beyond the abundance of their resources, like climate 
and natural enemies [28]. 

Saproxylics are strongly affected by the successional decay of 
wood, but they also contribute to determining the speed of this 
decay process. When saproxylic insects penetrate the suriiu c- of 
dead wood, they contribute to its physical degradation, improve 
aeration of the wood and facilitate the entry of decomposing fungi 
and bacteria [27,29]. This can increase the rate of decay of the 
wood [30,31]. Thus, the ability of saproxylic insects to respond 
numerically to the dead wood left by an insect outbreak may affect 
the speed with which the dead trees are decomposed, and, 
ultimately, the resilience of the whole forest ecosystem. Under- 
standing the numerical responses of saproxylic insects is, therefore, 
an important component of understanding the post-outbreak 
dynamics of forest ecosystems. However, most studies of how 
saproxylics respond to large-scale forest disturbances have dealt 
with fire and wind, whose effects are not necessarily direcdy 
comparable to those of insect outbreaks. Wind mainly produces 
lying dead trees, which tend to have high moisture content, while 
outbreaks leave mosdy standing dead trees, which are prone to 
desiccation. This difference is expected to cause contrasting 
successional trajectories in the saproxylic community, particularly 
in the early stages of succession [2 7] . Meanwhile, trees killed by 
fire favour a specialized guild of fire-adapted saproxylics [19]. As 
far as we are aware, the bark beetle outbreaks in the Bavarian 
Forest National Park in Germany represent the only case where 
the responses of saproxylic insects to insect outbreaks have been 
studied so far [24,32,33]. Thus, the knowledge of this subject is stiU 
very limited. 

In the present study, we examine the responses of saproxyUc 
beedes to a massive outbreak of the two geometrid moths 
(Lepidoptera: Geometridae) Epirrita autumnata (Bkh.) (autumnal 
moth) and Operophtem hrumata (L.) (winter moth) that occurred in 
the mountain birch [Betula pubescens ssp. czerepanovii Orlova) forest of 
the Varanger region in northern Norway (Fig. lA) during 2002— 
2009. Mortality of birch trees was high throughout this region and 
reached more than 95% locally [14,34]. Thereby, the availability 
of dead wood for saproxylic insects was elevated far above normal 
levels across a very large area. In order to assess the short-term 
responses of the saproxylic beede community in Varanger to this 
massive dead wood resource pulse, we used window traps to 
sample beetles in 2011 and 2012, both in severely damaged forest 
and undamaged forest at the cjutskirts of the outbreak area. Non- 
saproxylic beedes were also included in the study, in order to 
provide a contrast group that is not expected to benefit from 
increased availability of dead wood. By comparing beetle counts 
between damaged and undamaged forest, we address the following 
specific questions: 



1 . What is the magnitude of the short-term numerical responses 
of saproxylic beedes (both as individual species and as a 
community) to the superabundance of dead wood resources 
generated by the moth outbreak? 

2. Are the numerical responses of saproxylic beetles weak 
compared to the magnitude of dead wood resource input, 
suggesting that important limiting factors act to constrain 
beede responses? 

3. Have numerical responses occurred only in saproxylic beede 
species that are associated with dead wood in a particular stage 
of decay, suggesting decay-stage-preferences as a limiting factor 
for beetle responses? 

4. Has the outbreak caused changes in saproxylic beetle 
community structure (i.e. the relative abundances of different 
beede species), as would be expected if only a subset of the 
species in the community have responded numerically to the 
outbreak (cf question 3)? If so, what is the magnitude of this 
change in community structure? 

Materials and Methods 

Ethics statement 

No specific permissions were required to access any of the 
sampling locations or deploy any of the sampling methodologies 
described for this study. Three beede species that are included on 
the Norwegian red list [35] were incidentally caught (see below for 
description of the trapping method) in small numbers during the 
study, namely Atheta taxkeroide.s (lUCN category: Near threatened), 
Denticollis borealis (lUCN category: Vulnerable) and Enicmus 
lundbladi (lUCN category: Vulnerable). This has been reported 
to the Norwegain Biodiversity Information Centre (http://www. 
biodivcrsity.no/). Table SI in Appendix SI gives the number of 
individuals caught. 

Study region 

The Varanger region is located at approximately 70°N, 29°E, in 
the north-eastern part of Finnmark county, northern Norway 
(Fig. lA). Mean temperatures for January' and July are -12.2 and 
12.3°C, respectively, and annual precipitation is 400— 500 mm 
[Rustefjelbma meteorological station (70°23'55"N, 28°ir36"E), 
1961-2000 normal period]. The forest in the region is strongly 
dominated by mountain birch, although aspen (Populus tremula) and 
Scots pine {Pinus sylvestris) also occur patchily. Outbreaks by E. 
autumnata and 0. hrumata are the main natural disturbance factors 
in the north-Fennoscandian mountain birch ecosystem and occur 
approximately every 10 years [36-38], Outbreaks vary in 
intensity, but sometimes cause severe defoliation and high 
mortality rates of mountain birch across thousands of square 
km. The duration and severity of the last outbreak cycle was 
historically unprecedented, and it has been estimated that 
approximately 10 000 square km of birch forest in northern 
Fennoscandia was affected by severe defoliation at least once 
during the period 2002-2008 [34]. The Varanger region was 
particularly strongly affected due to consecutive severe outbreaks 
by E. autumnata (2002-2004) and 0. brumata (2005-2009). 

Study design and sampling 

In spring 2011, two transects were established in areas with 
slightiy different defoliation histories, namely Tana (70°03' N, 
27°45' E. Defoliated mainly during 2006-2007) and Kirkenes 
(69°46' N, 29°20' E. Defoliated mainly during 2007-2009) 
(Fig. IB). Both transects were roughly linear and consisted of 10 
sampling stations, interspaced by approximately two km. The first 
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Figure 1. Map of Finnmark County (A) with the study region in Varanger inset (B). Red areas represent birch forest which experienced 
severe defoliation during the most recent moth outbreak. Different shades indicate the main years of defoliation for different areas (see key). The 
sampling stations of the two transects are indicated by black dots. 
doi:1 0.1 371 /journal.pone.0099624.g001 



stations in each transect were located deep within forest which had 
suffered severe mortality due to the outbreak. Both transects then 
ran towards forest which was located at the outskirts of the 
outbreak area and had only experienced negligible defoliation 
(Fig. IB). Fig SI in Appendix SI shows photographs from the 
individual sampling stations. 

To assess the level of forest damage inflicted by the outbreak at 
the individual sampling stations, we established two 50 m survey 
lines at each station. The lines started at the centre-point of the 
station and ran towards the north and south, respectively. Every 
7 m along the lines, we selected the closest birch tree with at least 
one stem (standing or fallen) of more than 1.3 m in height. The 
damage of the three thickest stems (above a minimum height of 
1.3 m) within each selected tree was visually scored according to 
the criteria in table 1 . These registrations were only conducted in 
2011, since we had no reason to believe that the level of forest 
damage would change importantly between the two study-years 
(i.e. since the outbreak had ended in 2009). 

Inspection of damage scores for the individual sampling stations 
showed that the damage-level of the forest shifted abruptly along 
both transects (Fig. 2). Station 1-4 in Kirkenes and 1-5 in Tana 
had almost exclusively dead stems, and were therefore designated 
as "dead" stations. In contrast, at station 5-10 in Kirkenes and 7- 
10 in Tana, the great majority of stems were alive and either 
undamaged or with light loss of foliage. This is representative for 
the situation in a normal, healthy mountain birch forest. These 
stations were therefore designated as "living". Station 6 in Tana 
had a preponderance of dead stems, but also many living ones. 



The station was designated as dead, since its mean damage score 
fell closest to the dead stations. The simple dichotomy of dead and 
living stations was used to represent forest damage in all 
subsequent analyses. Levels of forest damage at the individual 
stations were representative for the surrounding landscape (Fig. IB. 
[34]). Thus, within a transect, dead and living stations were 
located within continuous belts of severely damaged and healthy 
forest, respectively. 

Sampling of beetles was conducted by means of large window 
traps, with two crossed 40x60 cm window panes [39]. AH traps 
were mounted between two adjacent birch stems, with the upper 
edge of the window panes approximately 1.5-2 m above the 
ground (Fig. SI in Appendix SI). We used three traps per 
sampling station (for a total of 60 traps across the two transects). 
One trap was located at the centre point of the station, while the 
other two were located approximately 50 m to the north and 
south, respectively. In both 2011 and 2012, the traps were 
mounted in early June and remained operative for approximately 
two months until early August (thereby covering the peak of insect 
activity in our high-latitude study system). 

All beetle individuals were identified to species level by an 
expert. Beetle species were thereafter classified according to their 
degree of association with dead wood, i.e. obligate saproxylic, 
facultative saproxylic or non-saproxylic, hereafter collectively 
referred to as dead wood association (DWA) groups. All species 
were also classified according to larval trophic guild, i.e. predators, 
fungivores, herbivores, wood-feeding, several (species belonging to 
more than one guild) or others (guilds that were rare in this study. 
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Table 1. Classification system used to score the damage of birch stems. 





Damage category 


Definition 


Damage score 


Living. Undamaged 


Stem with no apparent loss of foliage 


1 


Living. Lightly damaged 


Stem with >50% of foliage reaming compared to healthy stem 


2 


Living. Severely damaged 


Stem with <50% of foliage reaming compared to healthy stem 


3 


Dead 


Dead stem without foliage 


4 



doi:1 0.1 371 /journal.pone.0099624.t001 



i.e. scavengers, coprophages and detritus-feeders). Finally, sa- 
proxylic species were classified according to host-tree association 
i.e. associated with birch, not associated with birch and unknown 
host-tree association. The classification was based mainly on 
Dahlberg and Stokland [40] and Bohme [41], but additional 
information was obtained from a number of other sources. Table 
S 1 in Appendix S 1 gives a complete list of the sources used. 

Statistical analyses 

The goal of our analyses was to address the outlined questions 
about the numerical responses of saproxylic beetles in the study 
region to the changes in dead wood availability caused by the 
moth outbreak. The availability of dead wood was represented by 
our classification of forest damage, i.e. living forest (normal 
background levels of dead wood) and dead forest (greatiy elevated 
avaUabihty of dead wood owing the moth outbreak). Since there 
was an abrupt shift from living to dead forest along the transects, 
we also expected that beede communities would show rather sharp 
(i.e. non-gradual) spatial changes when moving from one forest 
type to the other. Therefore, we assumed that spatial changes in 
the beede community between living and dead forest were 
adequately modelled by a categorical living-dead contrast. The 
possibility of additional spatial patterns in the data was assessed by 
plotting station-specific beetie counts as well as topographical 
parameters along the transects. 

Forest damage (dead or living) was the focal predictor variable 
in all analyses. In addition, all models included the categorical 
predictors year (201 1 or 2012) and location (Kirkenes or Tana), to 



account for variation in beetie communities between study years 
and transects. Within each year, beetle counts were pooled across 
the entire trapping season and the three individual traps within 
each station, resulting in a single datapoint per station, per year. 
All analyses were conducted separately for the three DWA groups, 
owing to their different degrees of association with dead wood. For 
the two saproxyhc groups, the analyses only included species 
which are known to be associated with birch and which could 
therefore be expected to respond numerically to the dead birch 
wood left by the outbreak. 

We performed quantitatively oriented univariate analyses to 
assess the magnitude of the numerical responses of saproxylic 
beetles to the dead wood generated by the moth outbreak. First, 
we assessed whether the total abundances of the different DWA 
groups (pooled abundance across all species within each group) 
differed between dead and living forest. This was done by fitting 
log-linear models with the station-specific individual-counts of 
each DWA group as response variables, and forest damage, year, 
location and their two-way interactions as predictors. All models 
employed quasi-likelihood correction for overdispersion. We used 
the Akaike Information Criterion, corrected for small sample size 
and overdispersion (QAICJ, to select the best model for each 
DWA group. The minimal models included only the main effect of 
the focal predictor forest damage. 

We also investigated if the proportion that each DWA group 
constituted, out of the total number of individuals across all 
groups, differed between dead and living forest. This was done by 
fitting logistic regression models with the station-specific propor- 
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Figure 2. Station-specific damage scores of birch stems along the two sampling transects. The X axes represent the airline distances from 
the first sampling station in each transect. Large circles represent mean damage scores for each station, while small circles represent scores for 
individual stems. A small amount of noise has been added to the stem-specific scores to help distinguish individual datapoints. See table 1 for 
explanation of damage scores. 
doi:1 0.1 371/journal.pone.0099624.g002 
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tions of each DWA group as response variables. We employed the 
same predictor variables and model selection strategy for this 
analysis as for the analysis of the total counts of DWA groups 
(above). 

To assess numerical responses to dead wood availability in the 
most important saproxylic beetle species, we fitted log-linear 
models to the station-specific counts of all species which were 
judged to be sufficiently abundant for statistical modelling (total 
count >80 individuals). The predictor variables and model 
selection strategy for these species-specific models were the same 
as for the models fitted to the total counts of the DWA groups 
(above). 

Finally, we used multivariate methods to investigate how beetle 
community structure varied between living and dead forest. In 
order to assess overall community-responses to the outbreak, the 
multivariate models were run not only for the three DWA groups, 
but also collectively for all species in the data. All beetle counts 
were Hellinger transformed before multivariate analysis [42]. 

We performed redundancy analysis (RDA) (rda function in the 
R package vegan) to allow visualization and significance-testing of 
the effects of forest damage, year and location on beetle 
community structure. We assumed additive effects of the three 
predictors, since inclusion of interactions did not change the 
conclusions from the analysis. The significance of the marginal 
effects of the predictors (i.e. effects after accounting for all other 
predictors in the models) was assessed with permutational 
ANOVA (anova.cca function in the R package vegan) with a 
maximum of 9999 permutations. Site scores from the RDAs were 
displayed on biplots in order to investigate how the beede 
communities of the individual sampling stations were organized 
with respect to the predictor variables. 

Variance partitioning (varpart function in the R package 
vegan) was used to estimate how much of the variance in beetle 
community structure was explained by the focal predictor forest 
damage, as well as year and location. The results of the variance 
partitioning were visualised with venn diagrams, fitted with the 
venneuler function in the R package with the same name. 
Negative variance fractions were treated as zeros. 

Results 

Overall beetle community-structure 

A total of 5880 beede individuals, belonging to 207 different 
species, were caught during the study (Table 2. Table SI in 
Appendix SI provides a complete species list). ObUgate-, 
facultative- and non-saproxylic species accounted for 28.2%, 
12.5% and 58.1% of the individuals and 25.1%, 22.7% and 43.5% 
of the species, respectively. The degree of association with dead 
wood is unknown for 1.2% of the individuals and 8.7% of the 
species. Predators and fungivores were the most abundant and 
species-rich trophic guilds among obligate saproxylic species, while 
predators and herbivores were most abundant and species-rich 
among non-saproxylic species (Table 2). Predators also constituted 
the highest number of individuals and species among facultative 
saproxylics. Facultative saproxyhcs also included many individuals 
and species that were fungivores or belonged to several trophic 
guilds, as well as many individuals with unknown trophic guild 
(most of which belonged to the species Liotrichus affinis). Most of the 
obligate saproxylic individuals (85.5%) and species (78.8%) caught 
were associated with birch. For facultative saproxylics, most of the 
individuals (70.4%) and roughly half of the species (48.9%) were 
birch associated. 



Effects of dead wood availability on the total abundances 

and proportions of DWA groups 

Outputs from the log-linear models (see table S2 and S3 in 
Appendix S 1 for model selection and parameter estimates) for all 
three DWA groups were heavily influenced by large-scale spatial 
gradients in beetle counts along the transects (Fig. 3A-C). These 
gradients did not appear to have a consistent relationship with 
forest damage (i.e. their direction with respect to forest damage 
shifted between the two transects). Nevertheless, the results of the 
modelling were consistent with our expectation of a positive 
numerical response to increased dead wood availabiUty for 
obligate saproxylic species, but not for facultative saproxylic 
species. In particular, the predicted abundances of all three DWA 
groups were higher in dead than in living forest in Kirkenes. 
However, this effect was most pronounced for obligate saproxylics 
(Fig. 3A), which had roughly 3.5 times higher predicted 
abundance in dead forest in this transect, while the predicted 
abundances of facultative- and non-saproxyfics were only about 
1.5 and 2 times higher, respectively, in d("ad forest (Fig. 3B, C). In 
Tana, the patterns for facultative- and non-saproxylics showed a 
reversal with respect to forest damage, with predicted abundances 
that were roughly one third and half as high in dead as in living 
forest (Fig. 3B, C). Meanwhile, the predicted abundance of 
obligate saprox5'li('s was \'<'ry similar for the two forest types 
(Fig. 3A). These patterns \\<;re consistent across both study-years 
(i.e. model selection did not support damage x year interactions 
for any DWA group), although the overall abundances of all three 
groups were considerably higher in 2012. Thus, when the three 
DWA groups were contrasted, it would appear that obligate 
saproxylics were favoured by dead forest. This impression was 
supported by logistic modelling of the relative proportions of the 
DWA groups (see table S3 and S4 in Appendix SI for model 
selection and parameter estimates). The predicted proportion of 
obligate saproxylic individuals was about IQ% higher in dead than 
in living forest (odds ratio: 1.87. 95% CI: 1.22-2.87. Fig. 3D). 
Meanwhile, there was a roughly 4"/o decrease in the predicted 
proportion of non-saproxylic individuals in this forest type (odds 
ratio: 0.83. 95% CI: 0.59-1.16. Fig. 3E). These patterns were 
consistent both across locations and years (i.e. model selection did 
not support damage x year or damage x location interactions for 
these two groups). The predicted proportion of facultative 
saproxylic individuals was very similar in living and dead forest 
in 2011(odds ratio: 1.03. 95% CI: 0.51-2.07. Fig. 3F). However, 
the model selection supported the damage x year interaction for 
this DWA group, causing the predicted proportion of facultative 
saproxylics to be almost 10% lower in dead forest in 2012 (odds 
ratio: 0.39. 95% CI: 0.25-0.62. Fig. 3F). 

Effects of dead wood availability on individual beetle 

species 

Three birch associated obligate saproxylic, one birch associated 
facultative saproxylic and four non-saproxylic species were judged 
to be abundant enough for individual log-linear modelling (Fig. 4. 
See table S5 and S6 in Appendix SI for model selection and 
parameter estimates). These species collectively constituted 69.1% 
of all individuals in the data, GS.O'/n of all birch associated obligate 
saproxylic individuals, 81.2% of all birch associated facultative 
saproxylic individuals and 80.5% of all non-saproxylic individuals. 
The obligate saproxylic Alalthodes guttifer was also abundant enough 
for modelling, but not associated with birch. Nevertheless, we 
analysed this species for the sake of comparison with the other 
saproxylics. 
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Table 2. The total number of individuals and species (in brackets) caught during the study for each DWA group and larval trophic 
guild. 





Trophic guild 


Dead wood association (DWA) group 








Obligate sx. 


Facultative sx. 


Non-sx. 


Unknown 


Sum 


Predator 


597 [15] 


215 [25] 


1859 [44] 


29 [9] 


2700 [93] 


Fungivore 


820 [20] 


63 [14] 


65 [16] 


36 [6] 


984 [56] 


Herbivore 






1251 [13] 




1251 [13] 


Wood-feeding 


71 [5] 








71 [5] 


Several 


167 [11] 


33 [5] 


1 [1] 




201 [17] 


Other 






86 [9] 




86 [9] 


Unknown 


4 [1] 


422 [3] 


155 [7] 


6 [3] 


587 [14] 


Sum 


1659 [52] 


733 [47] 


3417 [90] 


71 [18] 


5880 [207] 



doi:l 0.1 371 /journal.pone.0099624.t002 



Two of the birch associated obligate saproxylic species, namely forest in both transects). Curiously, the non-saproxylic species 

Elateroides dermestoides (Fig. 4A) and Rahocerus foveolatus (Fig. 4B), Oxytelus laqueatus (Fig. 4F) showed a similar pattern (roughly 7.5 

showed evidence for a positive numerical response to dead wood times more abundant in dead forest in both transects). The 

availability, with higher predicted abundances in dead than living obligate saproxylic Podistra schoenherri (Fig. 4C) and the facultative 

forest in both Kirkenes and Tana [E. dermestoides: roughly 25 and 3 saproxylic L. affina (Fig. 4E) (both birch associated) did not appear 

times more abundant in dead forest in Kirkenes and Tana, to respond consistently to dead wood availability, with predicted 

respectively. R. foveolatus: roughly 2.5 times more abundant in dead abundances being higher in dead forest in Kirkenes and lower in 
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Figure 3. Variation in abundances and proportions of DWA groups along the two sampling transects. Large symbols: predicted 
abundances from log-linear models (A-C) and proportions from logistic models (D-F) in dead and living forest in Kirkenes and Tana in 201 1 and 2012. 
Error bars represent 95% confidence intervals. In cases where year effects were not supported by the model selection, only a single predicted value 
for each damage x location combination is presented, representing the overall estimate across 201 1 and 2012. Small symbols: original datapoints, i.e. 
the observed station-specific counts (A-C) and proportions (D-F) of obligate- facultative- and non-saproxylic beetle species in Kirkenes and Tana in 
2011 and 2012, plotted against the airline distance from the first station in each transect. 
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this forest type in Tana. A similar reversal of predicted abundance 
patterns between transects were evident for the not birch 
associated obligate saproxylic species M. guttifer (Fig. 4D) and the 
non-saproxylic species Anthophagus omalinus (Fig. 4G) and Eanus 
costalis (Fig. 4H). Finally, the non-saproxylic species Polydrusus 
Jiilvicomis (Fig. 41) was almost absent in dead forest but relatively 
abundant in Uving forest in both transects. 



Effects of dead wood availability on beetle community 
structure 

Forest damage explained none of the variance in community 
structure for facultative saproxylics (see below) and we therefore 
had no reason to expect that damage would yield significant effects 
in a RDA for this DWA group. Accordingly, we did not conduct 
RDA for facultative saproxylics. RDA site scores for living and 
dead forest tended to be separated in biplots for both obligate- and 
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non-saproxylic species, suggesting tliat these DWA groups showed 
consistent differences in community structure between the two 
forest types (Fig. 5). This separation was quite distinct for non- 
saproxylic species (Fig. 5B), but less clear for obligate saproxyhcs 
(Fig. 5A), with considerable overlap between living and dead 
stations for the latter group. Nevertheless, the effect of forest 
damage was highly significant in permutation-based ANOVAs for 
both DWA groups (/"= 0.005 in both cases). When all beetie 
species were considered collectively, the separation between living 
and dead stations became even clearer, with hardly any overlap 
between the two forest types (Fig. 5C). Also in this case, 
permutation-based ANOVA yielded a highly significant effect of 
forest damage {P= 0.005). The degree of separation between site 
scores for Kirkeiies and Tana and 2011 and 2012 was variable 
among the species groups, but all permutation-based ANOVA 
tests for location and year were significant at P= 0.030 or lower. 

Forest damage, location and year collectively explained 16.4% 
of the variance in beetle community structure for obligate 
saproxyhc species, 27.7% for non-saproxylic species and 20.1% 
for all species (Table 3, Fig. 6). For facultative saproxylics, damage 
failed to account for any variance, while location and year 
collectively explained 9.8%. In all cases, there was little overlap 
between the variance explained by the individual predictors, so 
that their total and unique contributions were largely similar. We 
therefore report only total contributions below. Damage was the 
most important predictor for obligate saproxylics, accounting for 
11.1% of the variance. Year was most important for facultative 
saproxylics, accounting for 7.4%. Year was also most important 
for non-saproxylics, accounting for 17.7%, while damage made a 
comparatively modest contribution of 7.6%. As expected from the 
results for the individual DWA groups, the relative contributions of 
year and damage became more even when all beetle species were 
considered collectively, with year explaining 10.2% and damage 
explaining 7.3%. Location made only a small contribution for all 
groups, explaining less than 4.0% of the variance in all cases. 

Discussion 

Insect outbreaks can rapidly generate more dead woody 
material than any other ecosystem process (with the possible 
exception of wind and fire), but littie is known of how saproxylic 
organisms respond to these disturbance events, despite the fact 
that saproxylics are intimately dependent on dead wood and play 
an important role in the ecosystem service of decomposition 
[43,44] . In this study we have contributed to filling this knowledge- 
gap by characterizing short-term changes in the abundance and 
community structure of saproxylic beetles following a severe moth 
outbreak in sub-arctic birch forest. 

We expected that the outbreak would generate a dead wood 
resource pulse and thereby trigger positive numerical responses in 
the saproxylic beetie community. Importantly, we also expected 
that such responses would be limited to a subset of the community, 
owing to the ability of given saproxyhc species to exploit only 
certain dead wood decay-stages, and that increases in responding 
species would be constrained by limitations to the amount of dead 
wood that can be colonized within a given amount of time (i.e. 
satiation effects). Our results appear to match these expectations. 
By contrasting the overall abundances of obligate-, facultative- and 
non-saproxylic species, we did indeed find evidence for increased 
abundance of obligate saproxyhcs in outbreak-affected forest. 
However, this increase corresponded only to a roughly 10% 
higher proportion of obligate saproxylic individuals in dead forest, 
while the increase in the dead wood resource may have been as 
much as hundredfold. 
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(corresponding to individuals sampling stations) from the RDAs. The 
letters on the symbols indicate the transect that each station belongs 
to. K= Kirkenes. T= Tana. Centroids for the levels of the three predictor 
variables are indicated by large letters on the plots. Forest damage: D = 
Dead. L= living. Location: Ki= Kirkenes. Ta= Tana. Yean 11=2011. 
12 = 2012.The two polygons in each plot surround all stations 
belonging to dead and living forest. Note that the site scores are 
shown as principal coordinates [65], implying that the distances 
between points are two-dimensional approximations of the Hellinger 
distances between them. The percentage of the total variation that is 
accounted for by the RDA axes is given in the axis labels. 
doi:10.1371/journal.pone.0099624.g005 

These fmding.s may .suggest that the saproxylic beetle commu- 
nity in the study region has a limited ability to respond to 
increased availability of dead wood. Our results for individual 
beetle species suggest species-specific decay-stage-preferences as 
one potential reason for this limitation. Most of the numerical 
response of obligate saproxylics was driven by the two species E. 
dermestoides and R. foveolatus. Both species exploit recentiy dead or 
dying trees and have therefore probably enjoyed increased 
resource-availability during the 3-5 years separating the conclu- 
sion of the outbreak and our study. Meanwhile, the third birch 
associated obligate saproxylic species that was analysed [P. 
schomheni) showed no response to the outbreak. This species 
utilizes dead wood in more advanced stages of decay and is 
thereby unlikely to have benefitted from the outbreak so far. Thus, 
our findings match the prediction that species-specific decay-stage- 



preferences will allow only a subset of the saproxylic beetle 
community to respond to an outbreak at any given time. 
Unfortunately, the knowledge of decay-stage-preferences for most 
species in our data is too limited to allow more rigorous evaluation 
of this hypothesis. 

Our study gives no direct indication of whether the numerical 
response of the saproxyhc beetie community is constrained by 
satiation effects. However, this seems likely to be the case. Within a 
given area, several decades can pass between moth outbreaks that 
cause substantial forest mortality. It is therefore unlikely that 
saproxylic beedes in the study region are specifically adapted to 
efficiendy exploit the extreme amounts of dead wood that 
outbreaks generate. Satiation effects could result e.g. from 
limitations to the amount of dead wood that adult beeties manage 
to colonize and reproduce in within a given time, including the 
number of eggs that they have available. Reproductive rates could 
be further constrained by the cold climate of the study region, 
which may prolong larval development of the decomposing insect 
community [45] . Such limitations on reproductive output may be 
decisive for the magnitude of numerical responses to the outbreak, 
since these responses are probably driven mainly by reproduction. 
Behavioural aggregation of beetles to the dead forest seems 
unlikely to be important, since the outbreak area is very large 
(Fig. 1 B) and a massive number of individuals would have to be 
drawn in from adjacent forest areas in order to raise densities 
appreciably. Moreover, the spatial extent of the outbreak area is 
probably far larger than the dispersal capacity of most beetle 
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Figure 6. Contributions of forest damage, location and year to explaining variance in beetle community structure. D= damage, L = 
location and Y= year. The amount of variance explained by each predictor is represented by the relative size of the circles in the venn diagrams. 
Overlapping areas indicate variance that is jointly explained by two of the predictors, while non-overlapping areas indicate variance that is uniquely 
explained by a single predictor. The residuals represent the proportion of the total variance in the beetle data that cannot be explained by the 
predictors. The values of the venneuler stress statistic [66] are below the 0.01 threshold (s.oi =0.056) in all cases, indicating that there is a good 
match (i.e. fit) between the venn diagrams and the actual variance components (Table 3). 
doi:1 0.1 371/journal.pone.0099624.g006 
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Table 3. Variance partitioning table showing how the predictor variables forest damage, location and year contribute to 
explaining the variance in beetle community structure for obligate-, facultative- and non-saproxylic (sx.) beetles, and for all beetle 
species collectively. 
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Location 
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0.038 


0.024 


0.035 


0.030 


Year 
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0.016 


0.074 
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0.102 


Damage + Location 
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0.149 


0.024 
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0.099 


Damage + Year 
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0.253 


0.175 


Location + Year 
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0.054 


0.098 


0.212 


0.132 


All 
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0.164 


0.098 


0.277 


0201 


Individual fractions 


Damage] Location + Year 
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0.111 
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0.063 


0.069 


Location | Damage + Year 
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0.038 


0.024 


0.022 


0.026 


Year | Damage + Location 


1 


0.016 


0.074 


0.177 


0.102 


Residuals 




0.836 


0.902 


0.723 


0.799 
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species in the data. Notably thougli, since the level of forest 
damage inflicted by the outbreak is spatially heterogeneous, 
aggregation on smaller spatial scale may affect local density 
patterns within the outbreak area. 

Climatic constraints on beetie diversity in our sub-arctic study 
region may also help explain the seemingly weak overall response 
of the saproxylic beetle community. In total, we encountered 52 
obligate- and 47 facultative saproxylic species. This is considerably 
less than what is typically reported from studies with comparable 
sampling efforts in more southerly deciduous forests [46-48] Thus, 
many species which may have responded to an outbreak further 
south may simply be absent from our system. Some support for 
this suggestion is provided by the results of Wermelinger et al. [21], 
who studied short-term responses of saproxyhc beetles to a large- 
scale windthrow event in Switzerland. Within 6 years of the 
windthrow, they observed much stronger increases in the 
abundance of saproxylic beetles than we did in the present study. 
This was largely due to species in the families Cerambycidae and 
Buprestidae and the subfamily Scolytinae (all of which contain 
birch associated species). These three groups were either absent 
(Buprestidae) or very rare (Cerambycidae: 1 species, 1 1 individ- 
uals. Scolytinae: 4 species, 16 individuals) in our data, presumably 
because our study region lies beyond the northern distributional 
limit of most of their species [25]. Our system may thereby lack 
many species that could potentially have responded numerically to 
a dead wood resource pulse. Another factor to consider is that our 
study region represents the northern distributional limit for most of 
the beetle species we encountered. Many species may thus be close 
to their tolerance limits for abiotic factors, like temperature, and 
could therefore be highly sensitive to stochastic year-to-year 
environmental variation [49] . This might also impair their ability 
to respond effectively to increases in resource availability. 

It should be noted that our interpretation of window trap 
catches in terms of beede numerical responses has some potential 
pitfalls. First, local sun exposure can potentially affect beetle 
activity levels, so that beetle counts are inflated in exposed (i.e. 
warmer) conditions. Sun exposure may also affect actual beetle 
abundance, since many saproxylic species have preferences for 
exposed or shaded conditions [50,51]. In our case, stations in dead 



forest (without foliage in the canopy) may have been more sun 
exposed than stations in living forest (with an intact canopy). 
However, total beetle counts were not consistently higher in dead 
forest (Fig. 3), suggesting that sun exposure did not inflate beetle 
counts in our case. Similarly, Kaila et al. [39] and Sverdrup- 
Thygeson and Ims [52] had similar total beede counts in exposed 
and shaded conditions. With respect to beetle habitat preferences, 
we are not aware that any of the most abundant species in our data 
(saproxylic or non-saproxylic) are considered to have strong 
preferences relating to sun exposure. Thus, we have no reason to 
expect that sun exposure was an important driver of the main 
abundance patterns in our study. Also, we should emphasize that 
even healthy mountain birch forest is quite open, with mostly small 
and well-spaced trees (This is well illustrated by Fig. SI in 
Appendix SI). Thus, the degree of sun exposure may not have 
differed radically between stations in living and dead forest. 
Second, saproxylic beetles may respond to variation in dead wood 
availability only on certain spatial scales, and window traps may 
fail to detect these responses if the volume of dead wood around 
the traps has not been quantified at the relevant scale [53-55]. 
However, in our case, the volume of dead wood was extremely 
high not only in the immediate vicinity of the traps, but also 
throughout the outbreak area as a whole, i.e. on a scale of 
thousands of hectares around the traps. Thus, in effect, our design 
captured elevated availabilities of dead wood on a wide range of 
spatial scales, ranging from locally around the traps and up to the 
landscape level. Thereby, we should have been able to detect 
beetie responses to the outbreak, regardless of whether beetles 
responded to dead wood availability on small or large spatial 
scales. Third, dispersal-driven "spillover" of saproxylic beetles 
from dead to living forest may have elevated beeties densities in 
living forest. This would invalidate our conclusion about weak 
numerical responses to the outbreak based on comparing beetie 
catches between the two forest types. However, such spillover 
would probably have been manifested as a gradual decline in 
saproxylic beetle counts when moving progressively further away 
from dead forest. This did not appear to be the case in our data. 
For beetle counts to be uniformly elevated throughout living forest 
there would have to be high rates of dispersal over distances of at 
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least 1 0 km (the approximate distance between the last living and 
dead stations in both transects). This seems improbable, although 
the dispersal capacity of the studied beedes is too poorly known to 
rule out the possibility either. Thus, in summary, although our 
results must be cautiously interpreted, we still consider it to be a 
reasonable conclusion that saproxylic beetles have so far respond- 
ed weakly to the dead wood left by the outbreak. 

If the ability of the saproxylic beede community to respond 
numerically to dead wood resource pulses is indeed limited, this 
could have important implications for the decomposition of trees 
killed by the outbreak. Saproxylic beetles contribute to wood 
decay by direct consumption and physical degradation of woody 
material, and by facilitating the entry and activity of decomposing 
fungi and bacteria [27,29-31,43,44]. However, if saproxylic 
beetles respond weakly to dead wood input, such decay processes 
may so far have played a minor role in decomposing the trees 
killed by the outbreak. Thereby, the decay of these trees may be 
retarded compared to dead wood in healthy mountain birch forest. 
This possibility must be explored in future studies, since we did not 
quantify the decay of dead stems along the transects. Nearly all 
stems were still standing in 2012 (pers. obs. by the authors), but 
this tells us little about decay rates, since birch stems have more 
than 80% probability of remaining standing throughout the first 
decade after their death (Makinen et al. [56] studying B. pubescens 
and B. penduh in southern and central Finland). Studies that 
contrast decay and insect-utilization of dead trunks between 
outbreak areas and healthy forest would be instructive, as would 
studies of fungal and other microbial decomposer communities. 
Penetration of bark by saproxylic invertebrates can be important 
for the colonization of dead trunks by bacteria and fungi [29], and 
weak numerical responses of saproxylic beedes to outbreaks may 
thereby retard colonization or alter community-composition in 
these microbial decomposers. Alternatively, the relative impor- 
tance of microbial decomposition may increase in the absence of 
strong numerical responses from saproxylic beetles. Studies 
addressing these questions would be of great interest for 
understanding post-outbreak decomposition dynamics. 

The indication of limitations in the short-term numerical 
response from the saproxyhc beetle community to increased dead 
wood amounts is also relevant in a management perspective. 
Several studies, especially in Northern Europe, have addressed the 
negative effects the decrease in background levels of dead wood in 
managed boreal forest have on the saproxylic beede communities 
[52, 57, .58]. In order to alleviate these effects, new management 
guidelines and restoration efforts are suggested to increase the 
dead wood amounts [53,59-61]. In this process, a firm under- 
standing of the responses of saproxyhc beedes across a wide span 
of spatial and temporal scales, and in different forest ecosystems, is 
vital. 

The separation of sites in living and dead forest in the RDA for 
obligate saproxylics indicates that the outbreak has induced 
changes in community structure (i.e. changes in the relative 
abundances of different species) among saproxylic beedes. This fits 
well with our expectation that only some saproxylic species should 
respond numerically to the outbreak, and the empirical observa- 
tion of species-specific response patterns. A more unexpected 
finding is the clearer separation (i.e. larger community changes) of 
dead and living sites for non-saproxylic species. This may suggest 
that there has been greater variation in species responses to the 
outbreak among non-saproxylics than obligate saproxyUcs. This 
seems plausible, since the non-saproxylic group contains more 
species and may therefore encompass a higher diversit)' of 
biological traits (e.g. microhabitat requirements and feeding 
specialisations). It is noteworthy in this respect that there were 



clear examples of both positive (i.e. 0. laqueatus) and negative (i.e. 
P. Julvicomis) responses to the outbreak among the individually 
analysed non-saproxylit: species, while this was not the case among 
saproxylic species. The mechanisms by which non-saproxylic 
species may have benefitted from the outbreak are unclear. With 
respect to 0. laqueatus, this species has been recorded in many types 
of decaying organic matter and also from sap runs [62]. Thus, we 
cannot rule out that this species has somehow benefitted from 
increased availability of dead wood, even though its substrate 
preferences seem to be too wide to justify classifying it as 
saproxylic. The finding of clear community changes in non- 
saproxylics means that our inferences about numerical responses 
to the outbreak, which were partiy based on contrasting the three 
DWA groups, must be treated with some caution. Clearly, non- 
saproxylics cannot be viewed as an unaffected control group with 
respect to the outbreak. Indeed, the pervasive impact of the 
outbreak on the whole ecosystem, including understory vegetation, 
herbivores [14,63] and abiotic conditions [64], probably leaves few 
biological taxa unaffected. However, this does not change our 
main conclusion about weak numerical responses in the saproxylic 
beede community. 

Although the RDA indicates that beetie community structure 
was affected by the outbreak, it should be emphasised that forest 
damage explained a small fraction of the overall variance in the 
beetle data (at most li.i"d, in the case of obligate saproxyhcs). 
Thus, although the outbreak represents a dramatic disturbance, 
which has changed both the biotic and abiotic environment, its 
effect on many beetle species may be o\erruled by other 
environmental variables. This conclusion is supported by the fact 
that counts of the DWA groups and many beetle species showed 
large-scale spatial gradients that appeared to be independent of 
forest damage. Notably though, these gradients may to some 
extent also represent variation in trapping efficiency or beetie 
activity, due to variation in e.g. wind speed or sun-exposure along 
the transects. Both transects have considerable heterogeneity in 
topographical variables, such as altitude and aspect (Fig. S2 in 
Appendix SI), and this may have affected both trapping 
efficiencies and beetie activity patterns. 

Conclusions 

Three to five years after an extensive moth outbreak, the 
saproxylic beetie community in the sub-arctic mountain birch 
forest appears to have mounted a weak but detectable numerical 
response to the increased availability of dead wood in the system. 
This response has mainly been driven by a few very abundant 
species, which utilize dead wood in early stages of decay. Although 
further monitoring is required to assess the fuU extent of beetie 
responses to the outbreak, our present results would suggest the 
possibility that saproxylic beetles increase in abundance to slowly 
to play a major role in decomposing th(- (;normous quantities of 
dead wood left by moth outbreaks in this system. Accordingly, the 
role of microbial decomposer communities in post-outbreak 
decomposition dynamics represents a prime topic for future 
research. 
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